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We present a combination of theoretical calculations and experiments for the low-lying vibrational 
excitations of H and D atoms adsorbed on the Pt(lll) surface. The vibrational band states are 
calculated based on the full three-dimensional adiabatic potential energy surface obtained from 
first-principles calculations. For coverages less than three quarters of a monolayer, the observed 
experimental high-resolution electron peaks at 31 and 68 meV are in excellent agreement with the 
theoretical transitions between selected bands. Our results convincingly demonstrate the need to go 
beyond the local harmonic oscillator picture to understand the dynamics of this system. 



Hydrogen on metal surfaces, and in particular on 
Pt(lll) has received considerable experimental [fj]-[7) and 
theoretical |j|-|l5| attention. Together with Pd, Pt is the 
most important material for heterogenous catalysis of hy- 
drogenation reactions. On the fundamental side, hydro- 
gen on metal surfaces provides a unique opportunity to 
observe the crossover from classical to quantum dynamics 
at relatively elevated temperatures [p~6] — p~9|| - 

For the vibrational excitations of H/Pt(lll), Baro 
et al. jl| applied high-resolution electron energy-loss 
spectroscopy (HREELS) and observed for H (D) at 
coverages up to 0.7 monolayer (ML) two peaks at 
68 (50) meV (assigned to the symmetric stretch mode 
SS) and 153 (112) meV (assigned to the asymmetric 
stretch mode AS). On the other hand, Richter and Ho 
H found three peaks at 67 (51) meV, 112 (84) meV and 
153 (108) meV in their HREELS measurements for a full 
ML of H (D). They assigned these to the AS mode, the 
SS mode, and the unresolved overtone and combination 
losses, respectively. To date, the disagreement between 
these observations is not resolved. In addition, recent 
measurement of the diffusion barrier for this system Q 
yielded a very low value of about 70 meV, considerably 
less than a previous value of about 200 meV. All these 
discrepancies emphasize our lack of understanding of this 
important system. 

On the theoretical side, an earlier first-principles (FP) 
study Q at a full ML coverage has shown that H adsorbs 
in the three-fold hollow fee site. More recent FP cal- 
culations pp|-p^| at lower coverages indicate that unlike 
for H/Ni(lll) and H/Pd(lll) the top site for 

H/Pt(lll) is also a local minimum only slightly higher in 
energy than the fee site. In addition, the energy barrier 
for the fcc-hcp-fcc path is only about 70 meV. In these 
calculations the main focus has been on adiabatic poten- 
tial energy surface (APES). The curvatures determined at 
the potential minima |p|,^0|JTT| could not account quanti- 



tatively for the observed vibrational excitation energies. 
Given the recent results of the shallow barrier between 
the fee and hep sites, it is clear that quantum effects are 
essential in understanding the dynamics of this system, 
and the vibrational spectroscopy need to be interpreted 
in terms of excitations between vibrational band states 
as proposed earlier [ p||l7j |. 

In this Letter we investigate the energetics and vibra- 
tional properties of H and D adatoms on the Pt(lll) sur- 
face. The experiments consist of HREELS investigations 
of H and D adsorption on Pt(lll) and its dependence on 
coverage. The experiments were performed in an ultra- 
high vacuum (UHV) apparatus with a base pressure of 
3 x 10~ n mbar as described elsewhere including sample 
preparation |24|]. The HREEL spectrometer was of the 
latest commercial Ibach design (Delta 0.5, SPECS, Ger- 
many). The H2 exposure was performed by backfilling 
the preparation chamber while the sample was cooled to 
85 K with liquid N2. The exposure was varied between 
0.1 and 500 L (1 L = 1.33 x 10" 6 mbarxs). The coverage 
was determined from thermal desorption (TD) assum- 
ing 9 = 1 at saturation. The TD spectra of Christmann 
et al. m were well reproduced. We focus on the spectrum 
for the low coverage interval 9 < 0.75 ML, within which 
the HREEL spectra are qualitatively the same |^E| . 

The spectra for H and D for 9 < 0.75 for specular ge- 
ometry are shown in Fig.|l]. It is well known that the loss 
intensity for vibrations of atomic H on transition metal 
surfaces is quite low since it is mainly caused by impact 
scattering. In spite of these intensity problems (note the 
magnifications in Fig. |l|) the spectra clearly indicate the 
two H derived losses at 31 and 68 meV. Note that these 
new observations are quite different from previous results 
obtained at a full ML [f|||]. In particular, the previously 
reported peaks at 112 meV and 153 meV are not observed 
at these lower coverages. In Fig. ||, we also show the 
results for off-specular geometry where the 31 meV peak 
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disappears while the 68 meV remains, indicating that the 
former is mostly of dipole character. For deuterium, there 
is a clear loss peak at 23.5 meV, corresponding to the 31 
meV peak for hydrogen with an isotope downshift. There 
is a very broad higher mode between 40 and 60 meV not 
very well resolved from the noise background. 

We now turn to the theoretical calculations for the vi- 
brational excitations. The FP calculations for 3D APES 
were done self-consistently with a parallel plane-wave 
code p6[ based on density functional theory with the 
GGA J27|. The Troullier-Martins pseudopotential Q 
was used for Pt and the Vanderbilt nonlocal ultrasoft 
pseudopotential [^9| for H. Electronic wave functions were 
expanded in a plane- wave basis with a kinetic energy cut- 
off at 35 Ry. For the irreducible Brillouin-zone integra- 
tion eight special k points |^| were used, together with 
a Fermi smearing of 0.2 eV. The calculated lattice con- 
stant for bulk Pt was a — 7.568 a.u. (experimentally 
a = 7.415 a.u. H). The supercell (2x2x4) has 16 
Pt atoms and one H atom, corresponding to a quarter 
monolayer coverage. The z axis is oriented towards the 
vacuum which was 14 A thick. The coordinates of the 
Pt atoms were relaxed before adding the H. Only the 
vertical coordinate of the H atom was relaxed during the 
calculations for each (at, y) position of the APES. We also 
checked the effects of the relaxation of Pt surface atoms. 
We found that the potential energy differences between 
the fee, bridge and hep sites changes less than 5 meV p2[ . 

We first present results for the 3D APES as shown in 
Fig. ||. The 2D minimum APES is shown in the inset. 
The top site is metastable and only 25 meV above the 
lowest energy adsorption fee site and is surrounded by 
large energy barriers. Along the fcc-hcp-fcc path, there 
is a shallow barrier of 78.9 meV. At all sites the H (D) is 
chemisorbed with an electronic charge transfer from the 
substrate toward the adatom, similar to that reported for 
H/Ni(lll) These results agrees well with other re- 

cent calculations ||l0|- |l3| , ^5l| . In particular, the approach 
of Kallen and Wahnstrom M| is closest to the present 
work and the resultant APES are very similar. How- 
ever, there are quantitative differences particularly con- 
cerning the curvature and anharmonicity at the fee sites. 
These will lead to substantial differences for the vibra- 
tional band excitations as detailed below. 

To study the vibrational modes of the H and D 
adatoms, we calculated the band states for the single 
particle Hamiltonian with the FP 3D APES. The Hamil- 
tonian is diagonalized in a basis set {|G, n)} given by a 
direct product between 595 in-plane plane waves (r\G) 
and 14 harmonic-oscillator states {z\^ n ) for the z coor- 
dinate. The lower branches n = 1 — 16 for H are shown 
schematically in Fig. |], along with the corresponding 
centers r*o of the probability distributions for the k = 
states, the band centers e n , and the widths As n . In addi- 



tion, many modes have mixed in-plane and vertical char- 
acter and cannot be described in conventional AS or SS 
terms. Also, the intervals between the band centers do 
not match with the local harmonic oscillator model based 
on the curvatures at the potential minimum. Finally, the 
higher bands have considerable bandwidth, reflecting the 
delocalized nature of the adsorbed H through tunnelling 
processes [ p|[l7| . 

The first band n = 1 localized at the fee site is very 
narrow (Aei < 0.1 meV). The next band n — 2 is local- 
ized at the hep site. Its center is at £2 = 23.6 meV above 
the center of the lowest band and it has a bandwidth of 
Ae2 = 1.3 meV. The difference £2 — £1 confirms that even 
with the zero point energies taken into account, adsorp- 
tion at low temperatures is confined to the fee sites. To 
analyze the vibrational excitation energies arising from 
transitions starting from the lowest band n = 1, we will 
focus only on those excited band states with finite ampli- 
tudes at the fee site. The first excitations are toward the 
composite band 304, which has a total bandwidth of 4.1 
meV. The two subbands are degenerate at k = 0, with E 
symmetry. The excitation to these bands, calculated as 
the difference between the bandcenters £394 — £1, occurs 
at 29.1 meV. Taking into account the bandwidths, this is 
in excellent agreement with the new mode at 31 meV ob- 
served in the present experiment. Moreover, this mode 
has mixed in-plane and vertical character with sizable 
dipole elements in the normal direction. This is confirmed 
by the fading of the 31 meV peak in the off-specular di- 
rection (see Fig.||). 

The next excitation is to the band n = 5, and to the 
7 © 8 pair of bands. The former is of A symmetry has a 
bandwith of 4.7 meV, while the latter has two subbands 
degenerate at k = (with E symmetry) and a total band- 
width of 11.2 meV. The excitation energies are 42 meV 
and 60 meV, respectively. These bands have mainly an 
in-plane character. We suggest that the observed peak 
centered at 68 meV could be a composite of these un- 
resolved transitions. Another point to note about the 
higher bands 7 ® 8 is that they have considerable weight 
at both the fee and hep sites. The states in these excited 
bands have high mobility throughout the surface. The 
excitation energy 60 meV needed to reach these bands 
is close to the activation energy observed in the quasi- 
elastic Helium Scattering experiment by Graham et al. 
p[. Thus, our band state results support the picture that 
the diffusion of H on Pt(lll) happens not through classi- 
cal hopping but rather via the activated tunnelling mech- 
anism as recently proposed |18| to account for the similar 
behavior observed for H/Ni(lll). The results give a di- 
rect confirmation of the protonic band picture introduced 
in @0. 

When D is substituted for H in the calculations, the 
lower bands have a one to one correspondence with the 
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H case with a downshift of frequency approximately by a 
factor of l/y/2 as expected. For deuterium, the hrst ex- 
citation energy is at 23.2 meV toward a composite band 
of width 0.2 meV, in excellent agreement with the ex- 
perimental observation for D as shown in Fig.|l|. There 
are transitions to higher bands at excitation energies be- 
tween 30 and 60 meV: an excitation to an A band at 
30.8 meV (bandwidth 0.8 meV), one to an E composite 
band at 46 meV (bandwidth 3.4 meV), and another to 
an E composite band at 59.8 meV (bandwidth 2.9 meV). 
These could account for the broad peak centered at 49 
meV observed in the experiment. 

There are many more excitations to the higher bands 
that we have not addressed here. Experimentally, they 
are not observed at low coverages (9 < 0.75 ML). Qualita- 
tively, one can understand this by noting that the higher 
bands are much more extended than the ground states 
and thus the matrix elements should be much smaller. 
As such, the resultant weak and broad peaks are hard to 
observe. The picture changes at full ML because then the 
tunnelling of the H atom would be suppressed by on-site 
exclusion and effectively the H atoms are localized again. 
Interestingly, at a full ML, the 31 meV peak disappears 
from the observations and two additional peaks appear 
at 110 meV and 153 meV Rj in the present experiment. 
These loss peaks agree nicely with the earlier results of 
Richter and Ho S. 

In summary, we have presented a combination of theo- 
retical calculations and experiments on the energetics and 
vibrational properties of H and D adatoms on Pt(lll). 
The experimental data for the loss peaks at low coverages 
(6 < 0.75 ML) differ significantly from previous studies 
at full H coverage. There is good agreement between 
the observed HREELS loss peaks and the protonic band 
structure calculations derived from a full 3D APES. This 
demonstrates that the vibrational spectroscopy of H ad- 
sorbed on metals and in particular for H/Pt(lll) cannot 
be understood simply in terms of localized AS and SS os- 
cillation modes at the adsorption sites, as often assumed 
in the literature. This picture is also crucial for a proper 
understanding of the recent diffusion data for H/Pt(lll) 
pL Finally, we emphasize that the dynamical behavior 
of this system at full ML is quite different from that at 
low coverages. There is also a possibility of change of 
adsorption sites. Further studies for the full ML case are 
in progress. 
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FIG. 2. Comparison of HREEL spectra for H on Pt(lll) 
for specular and off-specular geometries 



FIG. 1. HREEL spectra for H on Pt(lll) at 0=0.75 
ML(H) and 0.5 ML(D) for specular geometry. The spec- 
trum for D is in the inset, with the low-energy features dis- 
played on a smaller scale. The primary electron energy is 
Eo=2.0 eV. The exposure is given in Langmuir (L) (1 L 
« 1.33 x 10~ 6 mbarxs). 



FIG. 3. Vertical cuts in the 3D APES at high-symmetry 
sites. The topmost atomic plane is at z = 0, and the vacuum 
is at z > 0. The inset shows the 2D minimum APES along 
the path fcc-bridge-hcp-top-fcc. 
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FIG. 4. The lowest vibrational branches for H/Pt(lll). e n 
denotes the band center, Ae n the bandwidth, and ro the cen- 
ter of the BS with k = 0. 
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